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1  Introduction 

Plasmas  are  ionized  gas  mixtures,  either  magnetized  or  not,  that  have  many  practical 
applications:  lightning,  atmospheric  pressure  discharges  for  surface  treatment  and  flame 
stabilization,  Hall  thrusters  for  on-orbit  propulsion  tasks  on  communications  and  explo¬ 
ration  spacecraft,  laboratory  thermonuclear  fusion,  magnetic  reconnection  phenomenon  in 
astrophysics.  An  application  is  also  encountered  in  hypersonics  when  a  spacecraft  enters 
into  a  planetary  atmosphere  at  hypervelocity;  the  gas  temperature  and  pressure  strongly 
rise  through  a  shock  wave,  leading  to  excitation  of  the  particle  internal  energy  modes, 
dissociation  of  the  molecules,  and  ionization  of  the  particles  in  the  shock  layer  (Park, 
1990,  2008;  Tirsky,  1993;  Gnoffo,  1999;  Sarnia,  2000).  It  is  important  to  mention  that 
aerothermodynamic  heating  by  convection  and  radiation  dominates  the  design  of  the 
spacecraft  (Wright,  2008;  Prabhu,  2008).  Another  consequence  of  high  temperature  ef¬ 
fects  is  the  communications  blackout  experienced  during  a  part  of  the  entry  trajectory, 
where  is  it  impossible  to  transmit  radio  waves  either  to  or  from  the  vehicle.  Atmospheric 
entry  plasmas  are  reproduced  in  dedicated  wind-tunnels  such  as  plasmatrons,  arc-jet  fa¬ 
cilities,  and  shock-tubes  (Bogdanoff,  2008;  Bose,  2008;  Brun,  2008b;  Chazot,  2008;  Hollis, 
2008;  Panesi,  2008;  Surzhikov,  2008b). 

Depending  on  the  magnitude  of  the  ratio  of  the  reference  particle  mean  free  path  to 
the  system  characteristic  length  (Knudsen  number),  two  different  approaches  are  gener¬ 
ally  followed  to  describe  the  transport  of  mass,  momentum,  and  energy  in  a  plasma  (Bird, 
1994):  either  a  particle  approach  at  high  values  of  the  Knudsen  number  (see  Boyd,  2008, 
for  direct  simulation  Monte  Carlo  methods),  or  a  fluid  approach  at  low  values  (see  Can¬ 
dler,  2008;  Yee,  2008,  for  computational  fluid  dynamics  methods).  In  this  work,  we  study 
plasmas  that  can  be  described  by  means  of  a  fluid  approach.  In  this  case,  kinetic  theory 
can  be  used  to  obtain  the  governing  conservation  equations  and  transport  fluxes.  Hence, 
closure  of  the  problem  is  realized  at  the  microscopic  level  by  determining  from  experi¬ 
mental  measurements  either  the  potentials  of  interaction  between  the  gas  particles  or  the 
cross-sections,  allowing  for  the  transport  coefficients  to  be  computed.  We  do  not  consider 
the  transition  regime  between  the  continuum  and  rarefied  regimes  (see  Barth,  2008,  for 
Boltzmann  moment  methods). 

A  complete  model  of  plasmas  should  describe  the  following  physical  phenomena 

•  Thermal  nonequilibrium  of  the  translational  energy, 

•  Influence  of  the  electromagnetic  field, 

•  Occurrence  of  reactive  collisions, 

•  Excitation  of  internal  degrees  of  freedom. 

So  far,  no  such  unified  model  has  been  derived  by  means  of  kinetic  theory.  Besides,  a 
derivation  of  the  mathematical  structure  of  the  conservation  equations  also  appears  to  be 
crucial  in  the  design  of  the  associated  numerical  methods.  We  have  recently  investigated 
the  thermal  nonequilibrium  of  the  translational  energy  and  the  influence  of  the  electro¬ 
magnetic  field  (Graille  et  ah,  2009).  In  this  work,  we  summarize  the  main  results  for 
unmagnetized  and  weakly  magnetized  plasmas  relevant  to  hypersonic  applications.  Let 
us  now  describe  in  more  detail  how  these  issues  are  addressed  in  the  literature. 

First,  a  multiscale  analysis  is  essential  to  solve  the  Boltzmann  equation  governing  the 


RTO-EN-AVT-162 


2-3 


Kinetic  Theory  of  Plasmas 


ORGANIZATION 


velocity  distribution  functions.  We  recall  that  a  fluid  can  be  described  in  the  contin¬ 
uum  limit  provided  that  the  Knudsen  number  is  small.  In  the  case  of  plasmas,  a  thermal 
nonequilibrium  may  occur  between  the  velocity  distribution  functions  of  the  electrons  and 
heavy  particles  (atoms,  molecules,  and  ions),  given  the  strong  disparity  of  mass  between 
both  types  of  species.  The  square  root  of  the  ratio  of  the  electron  mass  to  a  character¬ 
istic  heavy-particle  mass  represents  an  additional  small  parameter  to  be  accounted  for 
in  the  derivation  of  an  asymptotic  solution  to  the  Boltzmann  equation.  The  literature 
abounds  with  expressions  for  the  scaling  for  the  perturbative  solution  method  (Devoto, 
1966;  Chmieleski  and  Ferziger,  1967;  Daybelge,  1970;  Kolesnikov,  1974;  Zhdanov,  2002). 
Petit  and  Darrozes  (1975)  have  suggested  that  the  only  sound  scaling  is  obtained  by  means 
of  a  dimensional  analysis  of  the  Boltzmann  equation.  Moreover,  they  have  deduced  that 
the  Knudsen  number  is  proportional  to  the  square  root  of  the  electron  heavy-particle  mass 
ratio.  Subsequently,  Degond  and  Lucquin-Desreux  (1996a, b)  have  established  a  formal 
theory  of  epochal  relaxation  based  on  such  a  scaling.  In  their  study,  the  mean  veloc¬ 
ity  of  the  electrons  was  shown  to  vanish  in  an  inertial  reference  frame.  Moreover,  the 
heavy-particle  diffusive  fluxes  were  scarcely  dealt  with  since  their  work  is  restricted  to  a 
single  type  of  heavy  particles,  and  thus  no  multicomponent  diffusion  was  to  be  found. 
In  such  a  simplified  context,  the  details  of  the  interaction  between  the  heavy  particles 
and  electrons  degenerate  and  the  positivity  of  the  entropy  production  is  straightforward. 
Here,  we  will  establish  a  theory  based  on  a  multiscale  analysis  for  multicomponent  plas¬ 
mas  (which  includes  the  single  heavy-particle  case)  where  the  mean  electron  velocity  is 
the  mean  heavy-particle  velocity  in  the  absence  of  external  forces.  As  an  alternative, 
Magin  and  Degrez  (2004b)  have  also  proposed  a  model  for  multicomponent  plasmas  in  a 
hydrodynamic  velocity  frame.  They  have  applied  a  multiscale  analysis  to  the  derivation  of 
the  multicomponent  transport  fluxes  and  coefficients.  However,  since  the  hydrodynamic 
velocity  is  used  to  define  the  reference  frame  instead  of  the  mean  heavy-particle  velocity, 
the  Chapman-Enskog  method  requires  additional  low  order  terms  in  the  integral  equa¬ 
tion  for  the  electron  perturbation  function  to  ensure  mass  conservation.  Finally,  we  also 
desire  that  the  development  of  thermal  equilibrium  models  shall  always  be  obtained  as  a 
particular  case  of  the  general  theory. 

Second,  the  magnetic  held  induces  anisotropic  transport  coefficients  when  the  electron 
collision  frequency  is  lower  than  the  electron  cyclotron  frequency  of  gyration  around  the 
magnetic  lines.  Braginskii  (1958)  has  studied  the  case  of  fully  ionized  plasmas  composed 
of  one  single  ion  species.  However,  the  scaling  used  in  this  study  does  not  comply  with  a 
dimensional  analysis  of  the  Boltzmann  equation.  Lucquin-Desreux  (1998,  2000)  has  inves¬ 
tigated  magnetized  plasmas  in  this  framework,  but  still  restricted  to  a  single  type  of  heavy 
particles.  Finally,  Giovangigli  and  Graille  (2003)  have  studied  the  Enskog  expansion  of 
magnetized  plasmas  and  obtained  macroscopic  equations  together  with  expressions  for 
transport  fluxes  and  coefficients.  Unfortunately,  they  did  not  account  for  the  difference  of 
mass  between  the  electrons  and  heavy  particles.  In  this  work,  we  will  deal  with  isotropic 
transport  coefficients,  typical  of  unmagnetized  plasmas,  such  as  atmospheric  entry  flows, 
and  weakly  magnetized  plasmas,  such  as  plasmatron  flows. 

Third,  plasmas  are  strongly  reactive  gas  mixtures.  The  kinetic  mechanism  comprises  nu¬ 
merous  reactions  (Capitelli  et  al.,  2000;  Bourdon,  2008;  Huo,  2008;  Laux,  2008;  Surzhikov, 
2008a):  dissociation  of  molecules  by  electron  and  heavy-particle  impact,  three  body  re¬ 
combination,  ionization  by  electron  and  heavy-particle  impact,  associative  ionization,  dis- 
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sociative  recombination,  radical  reactions,  charge  exchange. . .  Giovangigli  and  Massot 
(1998)  have  derived  a  formal  theory  of  chemically  reacting  flows  for  the  case  of  neutral 
gases.  Subsequently,  Giovangigli  and  Graille  (2003)  have  studied  the  case  of  ionized  gases. 
We  recall  that  their  scaling  does  not  take  into  account  the  mass  disparity  between  elec¬ 
trons  and  heavy  particles.  Besides,  in  chemical  equilibrium  situations,  a  long-standing 
theoretical  debate  in  the  literature  deals  with  nonuniqueness  of  the  two-temperature  Saha 
equation  for  quasi-equilibrium  composition.  Recently,  Giordano  and  Capitelli  (2001)  have 
emphasized  the  importance  of  the  physical  constraints  imposed  on  the  system  by  using  a 
thermodynamic  approach.  A  derivation  based  on  kinetic  theory  should  further  improve 
the  understanding  of  the  problem;  Choquet  et  al.  (2007)  have  already  studied  the  case  of 
ionization  reactions  by  electron  impact. 

Fourth,  molecules  rotate  and  vibrate,  and  moreover,  the  electronic  energy  levels  of  atoms 
and  molecules  can  be  excited.  Generally,  the  rotational  energy  mode  is  considered  to  be 
fully  excited  above  a  few  Kelvins.  In  a  plasma  environment,  the  vibrational  and  electronic 
energy  modes  are  also  significantly  excited.  The  detailed  treatment  of  the  internal  degrees 
of  freedom  is  however  beyond  the  scope  of  the  present  work  where  we  will  only  tackle  the 
translational  energy  in  the  context  of  thermal  nonequilibrium.  The  reader  is  thus  referred 
to  the  specialized  literature  for  the  treatment  of  the  internal  energy  (Brun,  2006,  2008a; 
Capitelli,  2008;  McCourt  et  ah,  1990;  Macheret,  2008;  Nagnibeda  and  Kustova,  2003; 
Schwenke,  2008). 

Fifth,  the  development  of  numerical  methods  to  solve  conservation  equations  relies  on 
the  identification  of  their  intrinsic  mathematical  structure.  For  instance,  the  system  of 
conservation  equations  of  mass,  momentum,  and  energy  is  known  to  be  nonconservative 
for  thermal  nonequilibrium  ionized  gases.  Therefore,  this  formulation  is  not  suitable  for 
numerical  approximations  of  discontinuous  solutions  (Liu  and  Vinokur,  1988;  Josyula  and 
Bailey,  2003).  Coquel  and  Marmignon  (1998)  have  derived  a  well-posed  conservative  for¬ 
mulation  based  on  a  phenomenological  approach.  Nevertheless,  their  derivation  is  not 
consistent  with  a  scaling  based  on  a  dimensional  analysis.  Our  kinetic  theory  derivation, 
based  on  first  principles,  naturally  allows  for  an  adequate  mathematical  structure  to  be 
obtained,  as  opposed  to  the  phenomenological  approach. 

In  this  work,  we  propose  to  derive  the  multicomponent  plasma  conservation  equations 
of  mass,  momentum,  and  energy,  as  well  as  the  expressions  for  the  associated  multicom¬ 
ponent  transport  fluxes  and  coefficients.  The  multicomponent  Navier-Stokes  regime  is 
reached  for  the  heavy  particles  and  is  coupled  to  first-order  drift-diffusion  equations  for 
the  electrons.  We  deal  here  with  first-order  equations  for  electrons,  thus  one  order  be¬ 
yond  the  expansion  commonly  investigated  in  the  literature.  The  derivation  relies  on 
kinetic  theory  and  is  based  on  the  ansatz  that  the  particles  of  the  plasma  are  inert  and 
only  possess  translational  degrees  of  freedom.  The  electromagnetic  held  influence  is  ac¬ 
counted  for.  In  Section  2,  we  introduce  the  Boltzmann  equation,  derive  the  H-theorem 
and  the  Maxwell  transfer  equations.  In  Section  3,  we  express  the  Boltzmann  equation  in 
the  heavy-particle  velocity  frame.  This  step  is  essential  to  establish  a  formalism  where 
the  electrons  follow  the  bulk  movement  of  the  plasma.  Then,  we  define  the  reference 
quantities  of  the  system  in  order  to  derive  the  scaling  of  the  Boltzmann  equation  from  a 
dimensional  analysis.  The  multiscale  aspect  occurs  in  both  the  streaming  operator  and 
collision  operator  of  the  Boltzmann  equation.  We  use  a  Chapman-Enskog  method  to 
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derive  macroscopic  conservation  equations.  The  system  is  examined  at  successive  orders 
of  approximation,  each  corresponding  to  a  physical  time  scale.  For  that  purpose,  scalar 
products  and  linearized  collision  operators  are  introduced.  The  global  expressions  for 
the  macroscopic  fluid  equations  are  then  provided  up  to  Navier-Stokes  equations  for  the 
heavy  particles  and  first-order  drift-diffusion  equations  for  the  electrons.  In  Section  4,  the 
multicomponent  transport  coefficients  are  calculated  in  terms  of  bracket  operators  whose 
mathematical  structure  allows  for  the  sign  of  the  transport  coefficients  to  be  determined, 
including  for  the  Kolesnikov  effect,  or  the  crossed  contributions  to  the  mass  and  energy 
transport  fiuxes  coupling  the  electrons  and  heavy  particles.  The  explicit  expressions  for 
the  transport  coefficients  by  means  of  a  Galerkin  spectral  method  (Chapman  and  Cowl¬ 
ing,  1939)  are  not  given  in  the  present  study.  Finally,  in  Section  5,  the  mass  conservation 
law,  the  first  and  second  laws  of  thermodynamics  are  proved  to  be  satisfied  by  deriving  a 
mass  conservation  equation,  a  total  energy  equation  and  an  entropy  equation.  Moreover, 
Onsager’s  reciprocal  relations  hold  between  the  transport  coefficients. 


2  Boltzmann  equation 

The  plasma  is  a  gas  mixture  composed  of  electrons,  denoted  by  the  index  e,  and  heavy 
particles  (atoms,  molecules,  and  ions),  denoted  by  the  set  of  indices  H;  the  full  mixture 
of  species  is  denoted  by  the  set  of  indices  S.  At  atmospheric  pressure  and  a  temperature 
of  0°  C,  a  perfect  gas  is  composed  of  3  x  1019  particles  in  a  cubic  centimeter.  Given 
the  enormous  number  of  particles  to  be  considered,  it  would  be  a  perfectly  hopeless  task 
to  attempt  to  describe  the  state  of  the  gas  by  specifying  the  so-called  microscopic  state, 

i.e.  the  position  x *  and  velocity  c*  of  every  individual  particle1,  and  we  must  resort  to 
statistics. 

2.1  Assumptions 

1.  The  plasma  is  described  by  the  kinetic  theory  of  gases  based  on  classical  mechanics, 
provided  that:  a)  The  mean  distance  between  the  gas  particles  l/(n°)1//3  is  larger 
than  the  thermal  de  Broglie  wavelength,  where  n°  is  a  reference  number  density 
(Hirschfelder  et  ah,  1954),  b)  The  square  of  electron  thermal  speed  14°  is  smaller 
than  the  square  of  the  speed  of  light. 

2.  Reactive  collisions  and  particle  internal  energy  are  not  accounted  for. 

3.  The  particle  interactions  are  modeled  as  binary  encounters  by  means  of  a  Boltzmann 
collision  operator,  provided  that:  a)  The  gas  is  sufficiently  dilute,  i.e.,  the  mean 
distance  between  the  gas  particles  l/(?r°)1/3  is  larger  than  the  particle  interaction 
distance  (u0)1//2,  where  a0  is  a  reference  differential  cross-section  common  to  all 
species,  b)  The  plasma  parameter,  quantity  proportional  to  the  number  of  electrons 
in  a  sphere  of  radius  equal  to  the  Debye  length,  is  supposed  to  be  large.  Hence, 
multiple  charged  particle  interactions  are  treated  as  equivalent  binary  collisions  by 

^Dimensional  quantities  are  denoted  by  the  superscript  *. 
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means  of  a  Coulomb  potential  screened  at  the  Debye  length  (Delcroix  and  Bers, 
1984;  Balescu,  1988). 

4.  The  plasma  is  composed  of  electrons  and  a  multicomponent  mixture  of  heavy  par¬ 
ticles.  The  ratio  of  the  electron  mass  m°e  to  a  characteristic  heavy-particle  mass  m°h 
is  such  that  the  nondimensional  number  e  =  {m°e/ is  small. 

5.  The  pseudo  Mach  number,  defined  as  a  reference  hydrodynamic  velocity  divided  by 
the  heavy-particle  thermal  speed,  Mh  =  v°/V®,  is  supposed  to  be  at  least  of  order 
one. 

6.  The  macroscopic  time  scale  t°  is  assumed  to  be  comparable  with  the  heavy-particle 
kinetic  time  scale  t°h  divided  by  e.  The  macroscopic  length  scale  is  based  on  a 
reference  convective  length  L°  =  v°t°. 

7.  The  reference  electrical  and  thermal  energies  of  the  system  are  of  the  same  order  of 
magnitude.  The  magnetic  field  intensity  is  such  that  the  electron  Hall  number,  or 
gyration  frequency  around  the  magnetic  line  multiplied  by  the  kinetic  time  scale,  is 
at  most  of  the  order  of  e. 

The  mean  free  path  1°  and  macroscopic  length  scale  L°  allow  for  the  Knudsen  number 
to  be  defined  Kn  =  l°/L°.  We  will  show  that  this  quantity  is  small,  provided  that 
assumptions  4-6  are  satisfied.  Therefore,  a  continuum  description  of  the  system  is  deemed 
to  be  possible. 


2.2  Inertial  reference  frame 


In  kinetic  theory,  the  plasma  particles  of  species  i  are  described  in  the  one-particle  phase 
space  (a?*,  c*)  by  f*,  their  velocity  distribution  function.  It  gives  the  probability  of  finding 
one  particle  of  species  i  at  position  x *  and  time  t*  with  a  velocity  c*.  In  other  words, 
^*da?*dc*  is  the  expected  number  of  i  species  particles  in  the  volume  element  da?*  located 
at  a?*,  whose  velocities  lie  in  the  velocity  element  dc*  about  velocity  c*  at  time  t*.  The 
gas  macroscopic  properties  at  that  position  are  then  obtained  as  average  microscopic 
properties  over  the  velocity.  The  partial  mass  density  p*,  the  hydrodynamic  velocity  v*, 
and  the  gas  thermal  energy  ev *  are  defined  by 


Pi 

p*v* 

p*ev *  +  |p*u*-u* 


I  f‘< 


dc* 


i  e  S, 

Ejesf 

£jeS///X*cTcid4 


(1) 


where  m*  is  the  mass  of  the  particle  i,  and  p*  =  p)  the  mixture  mass  density. 


Considering  assumptions  1-3,  the  temporal  evolution  of  the  velocity  distribution  func¬ 
tion  is  governed  by  the  Boltzmann  equation  (Chapman  and  Cowling,  1939;  Ferziger  and 
Kaper,  1972;  Cercignani  et  ah,  1994) 


a,.ff  +  c*-dx*f*  +  *;(«*  +  cp,B'ydc.f; 


£3*,(/*,p),  .es, 

3  6S 


(2) 
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where  symbol  E*  stands  for  the  electric  field,  B *,  the  magnetic  field,  and  q*,  the  charge 
of  the  particle  i.  The  left-hand-side  is  the  temporal  derivative  of  the  distribution  function 
c*,  t*)  in  the  phase  space 


ci  f*  cic  ■ 

Jr  =  ^*  +  c^*  +  ^.ac^*, 


accounting  for  the  acceleration  of  the  particles  due  to  the  Lorentz  force  q*(E*+c*AB*)/m* 
acting  on  the  charged  particles.  The  right-hand-side  is  the  collision  operator,  it  describes 
the  changes  of  the  distribution  function  due  to  interaction  with  all  the  particles  of  the 
mixture,  i.e.  collisions.  The  partial  collision  operator  of  particle  j  impacting  on  particle  i 
reads 


rv  (/.*■  //)  =  /  (/,*'//'  -  L'f;)  |c*  -  c*\a*jdujdCp  i,j  €  s. 


Quantities  after  collision  are  denoted  by  the  superscript ' .  The  particle  momentum  changes 
during  a  collision.  The  loss  term  is  proportional  to  the  number  of  particles  j  colliding  with 
particles  i,  thus  proportional  to  the  product  of  distributions  functions  if  we  assume 
a  lack  of  correlation  (molecular  chaos)  between  two  particles  that  are  about  to  collide. 
Considering  the  inverse  collision,  the  gain  term  is  proportional  to  It  is  necessary 

to  use  the  trick  of  the  inverse  collision  to  derive  the  expression  for  the  gain  term,  since  we 
can  not  assume  that  two  particles  that  have  just  collided  are  uncorrelated  (see  Chapman 
and  Cowling,  1939;  Ferziger  and  Kaper,  1972;  Cercignani  et  al.,  1994).  The  differential 
cross-section  cr*  is  defined  in  such  a  way  that  the  probable  number  of  collisions  per  unit 
volume  and  time,  for  which  the  direction  of  relative  velocity  after  collision  is  given  by 
io  =  (cf  —  c*')/\c*'  —  c*'|,  reads  —  c*|a*du;dc*dc*.  The  differential  cross-sections 

depend  on  the  relative  kinetic  energy  of  the  colliding  particles  and  the  cosine  of  the  angle 
between  the  unit  vectors  of  relative  velocities  lj  and  e  =  (c*  —  c*)/|c*  —  c*|.  The  partial 
collision  operator  of  particle  j  impacting  on  particle  i  is  obtained  by  integrating  over  all 
the  velocities  of  the  collision  partner  j  and  over  all  the  directions  of  the  relative  velocities 
after  collision.  We  notice  that  the  differential  cross  sections  are  symmetric  with  respect 
to  their  i,j  indices,  i.e.,  a*-  =  This  property  of  microreversibility  yields  the  sign  for 
the  entropy  production  rate,  implying  irreversibility  of  the  macroscopic  system. 


Boltzmann  gave  a  heuristic  derivation  of  his  equation  in  1872.  His  theory  was  attacked 
by  several  physicists  and  mathematicians  in  the  1890s,  because  it  appeared  to  produce 
paradoxical  results  (Ehrenfest  and  Ehrenfest,  1990).  However,  within  a  few  years  of 
Boltzmann’s  suicide  in  1906,  the  existence  of  atoms  had  been  definitively  established 
by  experiments  such  has  those  on  Brownian  motion.  Today,  Boltzmann’s  equation  is 
rigorously  derived  from  the  basic  laws  of  mechanics  in  the  so-called  Boltzmann-Grad 
limit  of  a  huge  number  of  particles  interacting  as  hard  spheres  of  very  small  diameter 
with  a  hnite  mean  free  path  value. 


2.3  H- Theorem 

Cercignani  (2006)  emphasizes  that  Boltzmann’s  motivation  when  deriving  his  equation 
seems  to  be  dealing  with  irreversibility  as  opposed  to  study  the  distribution  function  in 
time.  Boltzmann  has  also  introduced  a  quantity,  later  denoted  by  H,  defined  in  terms 
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of  the  velocity  distribution  functions.  He  then  established  that  as  a  consequence  of  his 
equation,  this  function  must  always  decrease  in  an  isolated  system  or,  at  most,  remain 
constant,  the  latter  case  occurring  only  if  a  state  of  statistical  equilibrium  prevails.  His 
result  is  usually  quoted  as  “//-Theorem”  and  indicates  that  H  must  be  proportional  to 
minus  the  entropy.  Based  on  this  result,  Boltzmann  showed  not  only  that  the  distribution 
introduced  by  Maxwell  is  a  steady  solution  of  his  equation,  but  that  no  other  solution  can 
be  found. 


We  consider  here  the  simplified  case  of  a  spatially  uniform  gas  in  the  absence  of  external 
forces  and  define  the  mathematical  entropy  H*  =  JAgS  /  f*  In  f*dc*.  By  taking  the  time 
derivative  of  this  quantity  and  using  the  Boltzmann  equation  for  this  simplified  case,  one 
obtains  after  some  algebra 


d/7* 

~dF 


a;)  M//T)  -  >"(/.*// 


|o-*-du>dc*dc*, 


as  shown  in  Appendix  A.  Since  the  function  (x  —  y)  (In  x  —  In  y)  is  non  negative,  d/P/dt*  < 
0  and  H*  is  a  non-increasing  function  of  time.  Based  on  energy  consideration,  it  can  be 
shown  that  H*  is  bounded  from  below  as  t*  —>  oo,  corresponding  to  a  state  in  which 
dH*/dt*  =  0  (Ferziger  and  Kaper,  1972).  Thus,  at  equilibrium,  the  integrand  vanishes, 
Si  Si  =  fi*fj *,  or  equivalently, 

In  fi  +  In  f*  =  In  f(*  +  In  f-*. 

This  property  is  typical  of  a  collisional  invariant  that  is  a  microscopic  quantity  globally 
conserved  during  a  collision  between  two  particles  i,j  G  S,  such  as  the  mass,  momentum, 
and  energy 


m*  =  mf, 


m *  =  m*', 


m*c *  +  m*c*  =  m*c *'  +  m*c *' 


3  3 


Introducing  the  species  vectors2 


^ + 4* 


«c-)ies’ 
(|m*c*-c*)  .gg, 


l  e  S, 

z/G  {1,2,3}, 


(3) 


where  Su  is  the  Kronecker  symbol  and  ?zs  denotes  the  number  of  species,  a  compact 
notation  is  found 

^  +  vf  =  Vf  +  V’f,  le{l,...,ns  +  4}. 

We  notice  that  the  mass  of  each  particle  is  conserved  during  the  collision.  Since  quantity 
In  f*  is  in  the  space  of  collisional  invariants,  it  can  be  expressed  as  a  linear  combination 
of  the  basis  vectors  i/j1* 


In  s:  =  m*A +  m*c*-A^*  + 

2Species  quantities  are  indicated  by  a  subscript.  When  the  subscript  is  dropped,  this  quantity  is 
a  species  vector.  Concerning  notation  of  spatial  quantities,  light-face  type  stands  for  spatial  scalars, 
bold-face  type  for  spatial  vectors  and  tensors. 
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where  the  coefficients  A\* ,  A 2*,  and  Af'  may  explicitely  depend  on  m*.  These  coefficients 
are  identified  by  means  of  eq.  (1).  For  that  purpose,  the  temperature  is  defined  in  kinetic 
theory  as 

|n*kBr«*  =  £  ///im*(c*  -  -  »*)  dc*.  (4) 

J6 S  J 


where  the  mixture  number  density  reads  n *  =  Ylj<=  s  ni>  wihli  the  partial  number  density 
71*  =  p*/m*.  The  superscript  v  indicates  that  temperature  is  measured  in  the  hydrody¬ 
namic  velocity  frame.  After  some  algebra,  we  obtain  the  Maxwell-Boltzmann  distribution 
function 


f?*  =  < 


m; 


2'k^T1 


3/2 


exp 


—m*(c*  —  v*)2 
2k  BTl 


i  e  S, 


(5) 


corresponding  to  the  velocity  distribution  for  a  gas  in  an  equilibrium  state. 


Boltzmann’s  //-Theorem  is  equivalent  to  the  second  law  of  thermodynamics  for  dilute 
gases.  It  is  important  to  mention  that  irreversibility  is  contained  in  the  collision  operator 
of  the  Boltzmann  equation  associated  with  the  molecular  chaos  assumption.  We  recall 
that  this  assumption,  needed  for  particles  that  are  about  to  collide,  is  rigorously  valid 
in  the  Boltzmann-Grad  limit.  It  can  be  shown  that  the  molecular  chaos  property,  if 
initially  present,  is  preserved  during  the  evolution  of  the  system.  Finally,  the  initial  chaos 
assumption  is  justified  as  follows.  Physically,  it  is  hard  to  prepare  an  initial  state  for 
which  the  molecular  chaos  assumption  does  not  hold.  The  physical  reason  is  that,  in 
general,  we  cannot  handle  the  single  particles  but  rather  act  on  the  gas  as  a  whole  at  a 
macroscopic  level,  usually  starting  from  an  equilibrium  state  (see  Cercignani  et  ah,  1994; 
Villani,  2002). 


2.4  Maxwell  transfer  equations 

The  transfer  equations  express  conservation  of  the  microscopic  properties  at  the  macro¬ 
scopic  level  of  the  flow.  They  are  obtained  by  multiplying  the  Boltzmann  equation  by  the 
collisional  invariants,  integrating  over  the  velocity,  and  summing  over  the  species  in  the 
mixture.  For  that  purpose,  let  us  introduce  the  scalar  product 

«r,or  =  E  /  dcv 

je  s  J 

for  families  £*  =  (£*)ieS  and  (*  =  s-3  It  is  important  to  mention  that  the  explicit 

expression  of  the  partial  collision  operator  3*j(.f*,  ./)*)  is  not  required  if  we  impose  that 
the  mass,  momentum  and  energy  are  conserved  in  the  mixture.  Actually,  Maxwell  had 
derived  transfer  equations  in  1867,  before  Botlzmann  proposed  an  expression  for  the 
collision  operator.  Projecting  eq.  (2)  on  the  collisional  invariants  given  in  eq.  (3),  one 
obtains 

+  {(c'-dx-r,  v>'*»  +  a A  (E‘ + cVB*)-»c*r,  v>'*»  =  o, 

m* 

3The  fully  contracted  product  in  space  is  denoted  by  symbol  “  0  ”,  for  instance,  for  two  scalar  a  and 
b ,  it  is  associated  with  a  product  ab ,  for  two  vectors  a  and  b,  a  scalar  product  a-b,  and  for  two  matrices 
A  and  B,  a  Frobenius  inner  product  A:B. 
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l  G  {1, . . .  ,  ns  +  4},  where  the  right-hand-side  vanishes.  The  conservation  equations  of 
mass,  momentum,  and  energy  are  derived  after  some  algebra 


dt*p*  +  dx*-(p*v*)  = -dx+-[J m*{c* -v*)f*dc*],  ie  S,  (6) 

dt*(p*v*)  +  dx*-(p*v*®v*)  =  -dx*-[J2  [ -  v*)f*dc*]  +  n*q*(E*  + 

jes  J 

dt*[p\ev *  +  §«*■«*)]  +  dx*-[p*v*(ev*  +  §«*■«*)]  =  -dx*-[J2  /  §Xci  ~  «2(c*  -  «*)# dc*] 

-X-[(^  J  m;(c*-i)>(c;-^)fdc;).^  +  (nY/+^;  J(c*-v*)f*dc*)-E\ 


(8) 


where  the  mixture  charge  is  given  by  n*q*  =  n*g*  and  symbol  ®  stands  for  the  tensor 
product.  The  expression  for  the  velocity  distribution  function  is  required  to  determine  the 
expression  of  the  transport  fluxes  and  close  the  system.  In  particular,  let  us  substitute 
the  Maxwcll-Boltzmann  distribution  eq.  (5)  for  f*  in  eqs.  (6)-(8) 


dt*P*  +  —  0,  *  G  S,  (9) 

dt*(p*v*)  +  dx*-(p*v*®v*)  =  -dx+-(p*I)  +  n*q*(E*  +  v*a  B*),  (10) 

dt*[p*(ev*+  lv*-v*)]  +  8x*-[p*v*(ev*  +  ±v*-v*)]  =  -dx*-(p*v*)  +  n*q*v*-E*,  (11) 

where  quantity  p *  =  n^k^T*  stands  for  the  mixture  pressure.  At  equilibrium,  the  fluxes 
of  mass  and  energy  vanish,  the  flux  of  momentum  is  the  hydrostatic  pressure  tensor, 
and  the  electric  current  is  only  due  to  convection;  the  transfer  equations  are  the  Euler 
equations.  In  order  to  describe  the  system  out  of  equilibrium  and  account  for  dissipative 
effects  (gradients  of  velocity,  temperature,  pressure,  and  species  concentrations)  in  the 
transport  fluxes,  a  perturbative  solution  method  can  be  used  to  linearize  the  Boltzmann 
equation  about  the  equilibrium  state  and  obtain  a  perturbed  solution  for  the  velocity 
distribution  function.  In  neutral  gases,  the  Knudsen  number  constitutes  the  small  per¬ 
turbation  parameter.  In  plasmas,  an  additional  small  perturbation  parameter  appears: 
the  square  root  of  the  electron  heavy-particle  mass  ratio.  This  parameter  is  essential 
to  describe  thermal  nonequilibrium  between  both  types  of  species,  associated  with  their 
quasi-equilibrium  states  at  the  electron  temperature  and  heavy-particle  temperature,  re¬ 
spectively.  A  multiscale  perturbative  method  for  plasmas  for  plasmas  is  proposed  in  the 
following  section. 


3  Multiscale  Chapman-Enskog  expansion 

In  about  the  same  year  (1916-1917),  Chapman  and  Enskog  independently  obtained  ap¬ 
proximate  solutions  of  the  Boltzmann  equation.  In  this  section,  we  adapt  the  Chapman- 
Enskog  perturbative  solution  method  to  the  multiscale  problem  of  the  Boltzmann  equation 
for  plasmas.  We  introduce  a  suitable  reference  frame  and  propose  a  scaling  of  this  equa¬ 
tion  based  on  a  dimensional  analysis.  A  translational  thermal  nonequilibrium  appears 
between  the  electrons  and  heavy  particles,  as  a  consequence  of  the  difference  of  mass 
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between  both  types  of  species.  Conservation  equations  are  derived  at  successive  orders  of 
approximation,  each  corresponding  to  a  physical  time  scale. 


3.1  Heavy- particle  velocity  frame 

Given  the  strong  disparity  of  mass  between  the  electrons  and  heavy  particles,  a  reference 
frame  linked  with  the  heavy  particles  appears  to  be  a  convenient  choice  for  plasmas.  This 
step  is  essential  to  easily  derive  a  formalism  where  the  electrons  follow  the  bulk  movement 
of  the  plasma.  We  define  the  mean  electron  velocity  and  mean  heavy-particle  velocity 


p>:  = 


[  m*ectf?dc*e,  p*hv*h  =  W 

J  je H  J 


m*c*//dc*, 


where  the  heavy-particle  mass  density  reads  p*h  =  Pj-  In  the  v*h  frame,  the  free 

electrons  interact  with  heavy  particles  whose  global  movement  is  frozen  in  space.  A 
similar  viewpoint  is  commonly  adopted  in  the  quantum  theory  of  molecules  when  the 
Born-Oppenheimer  approximation  is  used  to  study  the  motion  of  the  bound  electrons 
about  the  nuclei  (Born  and  Oppenhcimer,  1927).  Based  on  the  following  definition  of 
peculiar  velocities 

C*  =  c*-<,  i  G  S,  (12) 

the  heavy-particle  diffusion  flux  vanishes 


E 

16  H  ' 


m*c;^dc*  =  0. 


(13) 


Then,  the  Boltzmann  eq.  (2)  is  expressed  in  a  frame  moving  at  v*h  velocity  by  means  of 
this  change  of  variables 


+  (cr  +  O  -9**/*  +  Jr|E*  +  (C*  +  vl)  aB*]  -Sc;/.*  - 

-(Sq*/*®C*):8a!.<  =  y3*f(/i*,y),  us,  (14) 

16S 


where  the  material  derivative  reads  D/D t*  =  dt*  +  v^-dx*. 

3.2  Dimensional  analysis 

A  sound  scaling  of  the  Boltzmann  equation  is  deduced  from  a  dimensional  analysis  in¬ 
spired  by  Petit  and  Darrozes  (1975).  First,  reference  quantities  are  introduced  in  Table  1. 
The  characteristic  temperature,  number  density,  differential  cross-section,  mean  free  path, 
macroscopic  time  scale,  hydrodynamic  velocity,  macroscopic  length,  and  electric  and  mag¬ 
netic  fields  are  assumed  to  be  common  to  all  species.  Reference  dimensional  quantities 
are  denoted  by  the  superscript  °.  The  nondimensional  number 
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Common  to  all  species 

Temperature 

rpO 

Number  density 

n° 

Differential  cross-section 

a0 

Mean  free  path 

1° 

Macroscopic  time  scale 

t° 

Hydrodynamic  velocity 

v° 

Macroscopic  length 

L° 

Electric  held 

E° 

Magnetic  held 

B° 

Electrons 

Heavy  particles 

Mass 

m°e 

m°h 

Thermal  speed 

V° 

v  e 

vh° 

Kinetic  time  scale 

t° 

Le 

t°h 

Tabic  1:  Reference  quantities. 


quantifies  the  ratio  of  the  electron  mass  to  a  reference  heavy-particle  mass.  According  to 
assumption  4,  the  value  of  £  is  small.  Consequently,  electrons  exhibit  a  larger  thermal 
speed  than  that  of  heavy  particles 


K°  = 


/kRT° 


m\[ 


=  eV y 


(15) 


Moreover,  the  electron  and  heavy  particle  temperatures  may  be  distinct,  provided  that 
eq.  (15)  does  not  fail  to  describe  the  order  of  magnitude  of  the  thermal  speeds.  The 
differential  cross-sections  are  of  the  same  order  of  magnitude  a0.  Hence,  the  characteristic 
mean  free  path  1°  =  l/(n0cr°)  is  found  to  be  identical  for  all  species.  As  a  result,  the 
kinetic  time  scale,  or  relaxation  time  of  a  distribution  function  towards  its  respective 
quasi-equilibrium  state,  is  lower  for  electrons  than  for  heavy  particles 

,0  =  l_  ,0  =  l_  =  tj 

"  K0’  h  V«  e- 

Assumption  6  states  that  the  macroscopic  time  scale  reads 


t°  = 


(17) 


This  quantity  corresponds  to  the  average  time  during  which  electrons  and  heavy  particles 
exchange  their  energy  through  encounters.  In  addition,  the  macroscopic  temporal  and 
spatial  scales  are  linked  by  the  expression 


L° 


v°t°. 


where  the  hydrodynamic  velocity  is  determined  by  the  pseudo  Mach  number  Mj t 
Given  assumption  5,  the  Knudsen  number 


Kn 


Z° 


£ 

M~h 


V0/ VS- 


(18) 
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is  small,  due  to  our  choice  of  macroscopic  and  temporal  scales,  leading  to  a  continuum 
description  of  the  gas.  Finally,  following  assumption  7,  the  reference  electric  held  is  such 
that  q°E°L°  =  keT0.  The  intensity  of  the  magnetic  held  is  governed  by  the  Hall  numbers 
of  the  electrons  and  heavy  particles 


q°B° 


Pe  = 


rrut 


-t°  =  e1-^ 


a°B° 

ft,  =  q—t-L  =  eft. 

mV 


defined  as  the  Larmor  frequencies,  q°B°/me  for  the  electrons  and  q°B° /mQh  for  the  heavy 
particles,  multiplied  by  their  corresponding  kinetic  time  scale.  The  magnetic  held  is 
assumed  to  be  proportional  to  a  power  of  e  by  means  of  an  integer  b  <  0.  The  dimensional 
analysis  can  be  summarized  as  follows:  a)  Two  spatial  scales  were  introduced,  one  spatial 
scale  at  the  microscopic  level  and  one  spatial  scale  at  the  macroscopic  level,  they  are 
related  by  eq.  (18);  b)  Whereas  three  temporal  scales  were  defined  in  eq.  (16),  two  time 
scales  at  the  microscopic  level,  respectively  for  the  electrons  and  for  the  heavy  particles, 
and  one  time  scale  at  the  macroscopic  level,  given  in  eq.  (17),  common  to  all  species. 


Nondimensional  variables  are  based  on  the  reference  quantities.  They  are  denoted  by 
dropping  the  superscript  *.  In  particular,  one  has  the  following  expressions  for  the  particle 
velocities 

c*  =  He°ce,  c*  =  Vgd,  i  G  H. 

The  reference  hydrodynamic  velocity,  mean  electron  velocity,  and  mean  heavy-particle 
velocity  are  equal  to  v°.  We  investigate  the  system  at  the  macroscopic  time  t*  =  t°t 
and  macroscopic  length  x *  =  L°x.  The  reference  hydrodynamic  velocity,  mean  electron 
velocity,  and  mean  heavy-particle  velocity  are  equal  to  v°.  The  hydrodynamic  velocity  is 

(Ph  +  e2pe)v  =  phvh  +  e2peve, 

in  terms  of  nondimensional  variables,  whereas  the  mean  electron  and  heavy-particle  ve¬ 
locities  given  in  eq.  (3.1)  read 

peMhVe  =  -  f  Ce/edce,  phMhVh  =  ^  /  mjCjfjdCj. 

£  J  je h  J 

The  peculiar  velocities  defined  in  eq.  (12)  are  given  by  the  relations 
Ce  =  ce  -  eMhvh,  Q  =  d  -  Mhvh,  i  e  H. 

The  heavy-particle  diffusion  flux  vanishes,  as  shown  in  eq.  (13) 

[  "lJC’fjdCi  °-  (19) 

je H  ■' 

Thus,  the  Boltzmann  eq.  (14)  can  be  expressed  in  nondimensional  form  as 


dtfe  +  ^ia+eMbV^-d^+e-^q^a+sMhV^ABj-dcJe  +  i^qeE-eMn^-dcJe 

~  (dCJe®Cc):dxVh  =  ^[dee{fe,fe)  +  ^2dej(feJj)],  (20) 

j€H 
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dtfi  +  ^iQ  +  MhVHyd^+e^^Ka  +  MnVH^-dcJ  +  i^E-M^-dcJi 
~{dCi^i):da!vh  =  l[ldie{fiJe)  +  ^MJi)]>  i  e  H,  (21) 

je  h 

for  the  electrons  and  heavy  particles,  respectively.  Let  us  emphasize,  for  a  reader  familiar 
with  kinetic  theory,  that  eq.  (20)  for  the  electrons  exhibits  a  similar  scaling  as  that  of 
the  kinetic  equation  for  neutral  gases  in  the  low  Mach  number  limit  (yielding  parabolic 
macroscopic  equations),  whereas  the  scaling  of  eq.  (21)  for  the  heavy  particles  is  typical 
of  that  of  the  kinetic  equation  for  neutral  gases  in  the  compressible  gas  dynamics  regime 
(yielding  hyperbolic  macroscopic  equations).  Therefore,  the  coupled  system  of  kinetic 
eqs.  (20)  and  (21)  combines  the  usual  scalings  and  the  mathematical  structure  of  the 
resulting  system  of  macroscopic  equations  has  to  be  identified. 

The  multiscale  analysis  occurs  at  three  levels:  a)  In  the  kinetic  equations  (20)  and  (21); 
b)  In  the  collisional  invariants  and  thus  in  the  conservation  of  the  associated  macroscopic 
quantities;  c)  In  the  collision  operators.  Encounters  between  particles  of  the  same  type 
are  dealt  with  as  usual  (Chapman  and  Cowling,  1939),  whereas  the  collision  operators 
between  unlike  particles  (electron  heavy-particle  interactions)  depend  on  the  e  parameter 
via  their  relative  kinetic  energy  and  velocity,  and  the  vectors  cj  and  e.  The  scaling  of 
these  operators  is  investigated  in  Appendix  B. 

3.3  Enskog  expansion 

We  employ  an  Enskog  expansion  to  derive  an  approximate  solution  to  the  Boltzmann 
equations  by  expanding  the  species  distribution  functions  in  a  series  of  the  £  parameter 

fe  =  /e°(l  +  £4>e  +  4  +  ^K^), 

fi  =  i*°(i  +  £(k  +  £20j2)  +  CK^3),  f  €  h. 

Injecting  these  expressions  into  the  Boltzmann  eqs.  (20)  and  (20),  one  obtains 

£-l®e\fe°,  4  )  +  ^e°(/e°,  4  ,  <£)  +  (/e°,  4  ,  <t>l  <£) 

=  s-%2  +  e-1^-1  +  a°  +  eSl  +  0(£2),  (22) 

)  +  £®l(f?A)  =  £-1^  +  a°  +  4f!  +  0(e2),  i  e  H,  (23) 

given  in  more  detail  in  Appendix  C.  In  the  Chapman- Enskog  method,  the  plasma  is 
described  at  successive  orders  of  the  e  parameter  as  equivalent  to  as  many  time  scales. 
The  equations  derived  at  each  order  are  reviewed  in  Table  2.  In  the  remainder  of  Section 
3,  we  deal  with  the  macroscopic  equations,  in  Section  4,  with  the  expressions  for  the 
transport  fluxes. 

3.3.1  Thermalization 

We  solve  the  electron  Boltzmann  eq.  (22)  at  order  e~2 

ae~2  =  o, 
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Order 

Time 

Heavy  particles 

Electrons 

t° 

Le 

Expression  for  /e° 

Thermalization  (Te) 

e-1 

t° 

lh 

Expression  for  f®,  i  E  H 
Thermalization  (T/j) 

Equation  for  0e 

t° 

Equation  for  (j\ ,  i  E  H 
Euler  equations 

Equation  for  02 

Zero-order  drift-diffusion  equations 

£ 

t°/e 

Navier-Stokes  equations 

First-order  drift-diffusion  equations 

Tabic  2:  Chapman-Enskog  steps. 


corresponding  to  the  kinetic  time  scale  It  is  necessary  to  impose  mass  and  energy 
constraints 


(24) 


n,  =fffdC„ 

| n,T,  =  f  ff±C,-C,AC„ 

whereas  no  constraint  is  related  to  momentum.  The  electron  population  is  shown  to 
thermalize  in  the  heavy-particle  reference  frame  to  a  quasi-equilibrium  state  described  by 
a  Maxwell-Boltzmann  distribution  function  at  temperature  Te  (Graille  et  ah,  2009) 


3/2 

exp 


(25) 


Electrons  thermalize  in  any  velocity  frame,  however,  in  the  v ^  frame,  they  follow  the  bulk 
movement  associated  with  the  heavy  particles,  leading  to  a  physically  plausible  scenario. 
In  contrast,  heavy  particles  do  not  exhibit  any  ensemble  property  at  this  order. 


Then,  we  solve  the  heavy-particle  Boltzmann  eq.  (23)  at  order  e  1 

ar1  =  o, 

corresponding  to  the  kinetic  time  scale  t°h.  It  is  necessary  to  impose  mass,  momentum, 
and  energy  constraints 

rii  =  f  i  E  H, 

0  =  (26) 

\nhTh  =  ZwJ 

where  rih  is  the  heavy-particle  number  density.  The  heavy-particle  populations  are  shown 
to  thermalize  in  the  heavy-particlc  reference  frame  to  a  quasi-equilibrium  state  described 
by  a  Maxwell-Boltzmann  distribution  function  at  temperature  T \  (Graille  et  ah,  2009) 

f‘=ni{^k)  exp  (~whQ'c')  •  !£H’  (27) 
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The  quasi-equilibrium  states  are  described  by  means  of  distinct  temperatures  for  the 
electrons  and  heavy  particles.  This  is  not  contradiction  with  the  //-Theorem,  since  the 
equilibrium  state  has  not  been  reached  by  the  system. 


3.3.2  Linearized  collision  operators 

We  introduce  some  mathematical  tools  to  derive  the  conservation  equations  and  the  ex¬ 
pressions  of  the  transport  fluxes.  The  electron  linearized  collision  operator  fFe(0e)  and  the 
heavy-particle  linearized  collision  operator  3^  =  (3y)jeH  read 


3=e(4) 

Hth) 


—  __L 

/e°Je  ~  fe 


Z_^Ue3 

J'GH 


dee(fe%  ,  /e°)  +  0ee(/e°,  fe%  )  +  3°,(/e°4  ,  ff) 

i  e  H, 


— ^(3°  —  1 ,  fj°)  +  ,  f°4>j ) 

jSH 


where  /e°  is  given  by  eq.  (25),  ,/j°  by  eq.  (27),  and  (ph  =  (</>;);e h-  Then,  we  define  the 
electron  and  heavy-particle  scalar  products 


«fe,Ce»e=  /  fe©Ce  d Ce, 


ieH  J 


The  electron  collisional  invariants  are  introduced,  for  the  mass  and  energy,  as 

'i’l  =  i, 

^  =  ±Ce-Ce, 


(28) 


and  the  heavy-particle  collisional  invariants,  for  the  mass,  momentum,  and  energy,  as 

' '»!.  =  l  e  H, 

*■"+■'  =  (m,a„)i£H,  v  e  {1,2,3},  (29) 

>f+4  =  (| 

where  symbol  nH  denotes  the  number  of  heavy  particles  in  the  mixture.  The  linearized 
collision  operators  verify  the  following  properties 

«Te(0e),^»e  =  0,  /  e  {1,2}, 

«Wft),A»h  =  0,  l  e  nH  +  4}, 


for  all  functions  (f)e  and  (jfo . 


3.3.3  Electron  momentum  relation 

Macroscopic  equations,  such  as  the  transfer  eqs.  (6)-(8),  can  be  derived  by  means  of  a 
scalar  product.  The  projection  of  the  Boltzmann  eq.  (22)  at  order  e -1  on  the  collisional 
invariants  i(jle,  ^  *=  {1,  2},  is  trivial.  We  notice  that  momentum  is  not  included  in  the 
electron  collisional  invariants  given  in  eq.  (28),  as  opposed  to  full  set  of  collisional  invari¬ 
ants  given  in  eq.  (3).  This  is  due  to  the  structure  of  the  linearized  collisonal  operator 
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for  electrons  in  the  multiscale  analysis,  for  which  ((fFe(<^ ),  Ce))e  7^  0.  The  projection  of 
eq.  (22)  at  order  £_1  on  the  electron  momentum 

((®-\Ce))e  =  «ae-\  Ce»e, 


yields  a  macroscopic  relation  for  the  zero-order  momentum  transfered  from  electrons  to 
heavy  particles 


jSH 


1  «  Mep 

- CDWp - hj, 

Mh  Mk  ’ 


(30) 


where  the  electron  pressure  reads  pe.  This  original  property  of  the  Chapman-Enskog  ex¬ 
pansion  at  this  order  is  associated  with  the  absence  of  a  momentum  constraint  in  eq.  (26). 
Finally,  the  net  zero-order  momentum  exchanged  between  electrons  and  heavy  particles 
vanishes,  i.e. 


«3LU°,/jv o.V'jy"}}, 

for  v  G  {1,  2, 3}. 


jeH 


(/ev=y0).cs„))«  =  o, 


(31) 


3.3.4  Zero-order  electron  drift-diffusion  and  Euler  equations 

We  derive  zero-order  electron  drift-diffusion  equations  and  Euler  equations  based  on  the 
Boltzmann  equations  at  order  e°  corresponding  to  the  macroscopic  time  scale  t°.  We 
project  eqs.  (22)-(23)  on  the  collisional  invariants, 

«9?,#»e  =  «MX  ^  e  {1,2}, 

«,*»,  =  ((3Ulh))h,  ZG{l,...,nH  +  4}. 

After  some  algebra  and  using  eqs.  (30)-(31),  one  obtains 


dtpe  ”t“  9 x' 


Pe(Vh  +  jfcK) 


0, 


ft(PeCe)  T  9x'  h)  Pc^x'^h  T  ^^01 

dtpi  +  dx-(piVh)  =  0,  ie  H, 

dt(phVh)  +  dx-(phvh®vh  +  -jrppt)  =  jjpnqE  +  Sb0I0AB, 

1h  1h 

dt(phZh)  +  dxiphZhVh)  -  -Vhdx-Vh  +  A E%, 


(32) 

(33) 

(34) 

(35) 

(36) 


where  the  heavy-particle  pressure  is  ph  =  nhTh,  and  the  mixture  pressure,  p  =  pe  +  ph- 
Dissipative  effects  already  appear  at  this  time  scale,  we  introduce  the  electron  diffusive 
velocity,  the  electron  conduction  current  density,  the  total  current  density 


ve  =  -  f  CJe%  dCe,  Je  =  neqeVe,  Io  =  nqvh  +  d-  Je,  (37) 

ne  J  Mh 

and  the  electron  heat  flux 

qc  =  I  lce.cecJe%ace.  (38) 

It  is  important  to  mention  that  the  electrons  diffuse  in  the  heavy-particle  velocity  frame, 
following  the  bulk  movement,  even  though  they  thermalize  in  any  reference  frame.  Finally, 
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the  energy  transferee!  from  heavy  particles  to  electrons  at  order  zero  follows  a  Landau- 
Teller  relaxation  term 


A  E%=* 


1  ne(Te  -  Th 


T 


E 


2  n 


3 


o  ^7'ej 

inem~ 
je  H  e 


where  r  is  the  relaxation  time.  The  expression  for  the  collision  frequencies  uie  is  given 
in  Appendix  B.  The  net  zero-order  momentum  exchanged  between  electrons  and  heavy 
particles  vanishes 

A  El  +  A  E°h  =  0. 


3.3.5  First-order  electron  drift-diffusion  and  Navier-Stokes  equations 

We  derive  first-order  electron  drift-diffusion  equations  and  Navier-Stokes  equations  based 
on  the  Boltzmann  equations  at  order  e  corresponding  to  the  macroscopic  time  scale  t°/e. 
We  project  eqs.  (22)-(23)  on  the  collisional  invariants, 


«^e))e  +  ^e^e))e  =  <<«X  +  *<<&  #>>.>  «  ^  {1,2}, 


^h^lh))h  +  £i^h^h))h  =  +  le{l,...,nn  +  A}. 

After  some  algebra,  one  obtains 

dtpe  +  dx-  pe(vh  +  j^(Ve  +  eVe2))  =0, 

dt(peee)  +  dx •  ( peeevh )  =  -pedx-vh  -  A- dx ■  (qe  +  eq2c)  +  ^  (Je  +  £j‘e  )  E 
+  SbozMfrJe-VhAB  +  A  E®  +  eAFg, 


dtPi  +  d x-  Pi(vh+  ^-Vi) 


=  0,  i  e  H, 


(39) 

(40) 

(41) 

(42) 


dt(phVh)  +  dx-{phvh®vh  +  Ajpl)  =  ~^dx-nh  +  j^nqE  +  5wI0aB , 

±h  ±h  ±h 

dt(pheh )  +  dx-(phehvh)  =  ~(ph I  +  ellh)\dxvh  -  j^dx-qh  +  ^  Jh  E  +  A E%  +  eA Elh. 


Mh' 


(43) 


Like  in  the  case  of  neutral  gases,  dissipative  effects  for  the  heavy  particles  appear  at 
this  time  scale.  We  introduce  the  heavy-particle  diffusion  velocities,  the  heavy-particle 
conduction  current  density 


Vi  =  ~  I  dCi,  Jh  =  1  e  H> 

Ui  J  JGH 

the  heavy-particle  viscous  tensor, 


and  the  heavy-particle  heat  flux 


(44) 


(45) 


(46) 
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For  a  multicomponent  mixture  of  heavy  particles,  the  model  for  the  electron  transport 
fuxes  is  extented  one  order  further  than  in  the  literature.  We  introduce  the  second-order 
electron  diffusive  velocity,  the  second-order  electron  conduction  current  density 


Ve2  =  -  Ce/eVe  dCe,  J e2  =  neqeVc2 , 


nP 


and  the  second-order  electron  heat  flux 


<  =  /  -C.-C.C.j^dC.. 


(47) 


(48) 


The  energy  transfered  from  heavy  particles  to  electrons  at  order  one 

A  El  =  Y,n’VfFi”  («) 

je  h 

is  expressed  by  means  of  the  average  force  of  an  electron  acting  on  a  heavy  particle  i 


Fie  =  I  QS’dT.r)  I7el7e  fe( 7e) A  <7e)  d7e.  *  £  H, 


(50) 


where  the  expression  for  the  average  momentum  cross-section  Q^  '  is  given  in  Appendix  B. 
The  net  first-order  momentum  exchanged  between  electrons  and  heavy  particles  vanishes 

AEl  +  A  El  =  0. 

The  second-ordre  electron  transport  fluxes  given  in  eqs.  (47)- (48)  and  the  first-order  energy 
echange  term  given  in  eq.  (49)  are  essential  to  derive  a  total  energy  equation  and  an  entropy 
equation  that  satisfy  respectively  the  first  and  second  laws  of  thermodynamics. 


4  Transport  fluxes  and  coefficients 

In  this  section,  we  derive  the  perturbation  functions  yielding  the  expressions  for  the  trans¬ 
port  fluxes  that  appear  in  the  macroscopic  equations.  We  give  a  complete  description  of 
the  Kolesnikov  effect,  or  crossed  contributions  to  the  mass  and  energy  transport  fluxes 
for  multicomponent  plasmas  coupling  the  electrons  and  heavy  particles. 


4.1  First-order  electron  perturbation  function 


The  first-order  electron  perturbation  function  </>e  is  a  solution  to  the  electron  Boltzmann 
eq.  (22)  at  order  e_1,  he., 


2e(&)  = 


~Pe 


4'De  cL  -  d 


Ai 


where  4r(?e  and  are  functions  of  Cc  (Graille  et  ah,  2009).  The  electron  diffusion 
driving  force  is  defined  as  de  =  (dxpe  —  neqeE) / pe.  Uniqueness  of  the  solution  is  imposed 
by  means  of  the  constraints 

(</eVe,^»e  =  0,  Z  G  {1,2}, 
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i.e.,  the  perturbation  function  does  not  contribute  to  the  macroscopic  electron  mass  and 
energy.  By  linearity  and  isotropy  of  the  linearized  collision  operator  the  electron 
perturbation  function  is  expressed  in  terms  of  the  electron  driving  force  and  temperature 
gradient 

=  -Pe^eC-dc  -  0^-0®  (^r),  (51) 

where  the  vectorial  functions  (f)^e  and  <f>^e  are  the  solutions  to  the  equations 


x  e  i 


under  the  scalar  constraints 


with  n  G  {Dc,  Ag}. 


«/,°0^»e 


o,  l  e  {1,2}, 


The  transport  fluxes  are  obtained  by  injecting  the  expression  for  given  in  eq.  (51) 
into  the  eqs.  (37)  and  (38)  and  introducing  the  electron  bracket  operator 

Be,  Cl  =  ((fX.UQ)),- 

We  obtain  the  electron  diffusion  velocity  and  heat  flux 

Ve  =  —Dede  -  9edx\nTe,  qe  =  -\'edxTe  -  pe9ede  +  peheVe ,  (52) 


where  the  electron  diffusion  coefficients,  thermal  diffusion  coefficients,  and  partial  thermal 
conductivity  are  given  by 


D ,  =  0ea]„  ».  =  -p4I<Aea,  <A]„  a;  = 


The  first  term  of  the  diffusion  velocity  in  eq.  (52)  yields  diffusion  effects  due  to  the  partial 
pressure  gradient  and  electric  field.  The  second  term  represent  diffusion  arising  from  the 
electron  temperature  gradient  and  is  termed  the  Soret  effect.  The  first  term  of  the  heat 
flux  represents  Fourrier’s  law,  the  second  term  corresponds  to  the  Dufour  effect,  that  is, 
heat  diffusion  due  to  the  electron  partial  pressure  gradient  and  electric  field,  which  is  the 
symmetric  of  the  Soret  effect.  Finally,  the  third  term  is  the  transfer  of  electron  energy 
due  to  diffusion.  Alternative  forms  of  the  electron  diffusion  velocity  and  heat  flux  are  also 
introduced 


Ve  =  - De(de  +  xe9JnTe),  qc  =  -XcdxTe  +  peXeVe  +  PeheVe,  (53) 

where  the  thermal  diffusion  ratio  ye  is  defined  by  the  relation  9e  =  Deye,  and  the  thermal 
conductivity  by  Ae  =  A'e  —  ne\e9e. 

4.2  First-order  heavy-particle  perturbation  function 

The  first-order  heavy-particle  perturbation  function  =  (<A);£h  a  solution  to  the 
heavy-particle  Boltzmann  eq.  (23)  at  order  e°,  i.e., 

U4>h)  =  £HW)  +  3?] 

=  -  p,,  W  vPiJj  -  (v)  -  Tf  (T,  -  Th),  i  e  H, 

ie  H 
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where  /i  G  {qh,  X'h,  0},  are  functions  of  C)  (Graille  et  ah,  2009).  A  linearly 

independent  family  of  diffusion  driving  forces  is  also  introduced 


d , 


if1 - 


nfli 

Ph 


E 


rijMh  pi 
Ph  ie’ 


i  g  H, 


(54) 


where  quantity  pt  =  riiTh  stands  for  the  partial  pressure  of  species  i  6  H.  Uniqueness  of 
the  solution  is  imposed  by  means  of  the  constraints 

((fUh,  fth))h  =  0,  l  G  {1, . . . ,  nH+4}, 


he.,  the  perturbation  function  does  not  contribute  to  the  macroscopic  heavy-particle  mass, 
momentum,  and  energy.  By  linearity  and  isotropy  of  the  linearized  collision  operator 
T),  the  heavy-particle  perturbation  function  is  expressed  in  terms  of  the  heavy-particle 
velocity  gradient,  diffusion  driving  force,  temperature  gradient,  and  the  difference  between 
the  electron  temperature  and  heavy-particle  temperature 

<k  =  -4>T:d*vh  -  Ph  4>?-dj  -  <f>ih-9x  -  (j)f  (Te  -  Th),  i  G  H,  (55) 

je  H 

where  the  families  of  tensorial  functions  =  (0^)ieH,  of  vectorial  functions  (f>^3  =  (0fJ)jGH, 
j  G  H,  and  <phh  =  h,  and  of  scalar  functions  0®  =  (0f)igH  are  the  solutions  to  the 

equations 

H<K)  =  *?,  *  e  h, 

under  the  scalar  constraints 


0,  l  G  {1, . . . ,  nH  +  4}, 


with  \i  G  {rjh,  {Dj)je H,  A),,  0}. 


The  transport  fluxes  are  obtained  by  injecting  the  expression  for  (f\  given  in  eq.  (55) 
into  the  eqs.  (44)- (46)  and  introducing  the  heavy-particle  bracket  operator 


[Zh,Ch]h  =  ((fith,3:h(Ch)))h- 

We  obtain  the  diffusion  velocity  of  species  i  G  H 


v  =  -  E  d‘A  - 

jeH 


(56) 


where  the  diffusion  coefficients  and  thermal  diffusion  coefficients  are  given  by 

=  >,je  H,  =  ,eH. 


The  diffusion  matrix  DV]  is  symmetric.  The  first  term  of  the  diffusion  velocity  in  eq.  (56) 
yields  diffusion  effects  due  to  the  partial  pressure  gradient,  electric  field,  and  average  elec¬ 
tron  force.  The  second  term  represent  diffusion  arising  from  the  heavy-partcile  tempera¬ 
ture  gradient  (Soret  effect).  It  is  interesting  to  notice  that,  before  the  work  of  Chapman 
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and  Enskog,  thermal  diffusion  in  the  gas  phase  had  been  unknown  theoretically  and  unob¬ 
served  experimentally.  Then,  we  introduce  the  tensor  S  =  [dxVh  +  (dxVh)T]  —  \dx-Vh  I. 
The  heavy-particle  viscous  tensor  reads 

nh  =  —r)hs,  (57) 

where  the  shear  viscosity  is  given  by  rjh  =  Th{4>nhh,  ]/(/10.  The  heavy-particle  heat  flux 
reads 

q,  =  -A  'hd«Th  »A  +  E  PihiVi  ■  (5S) 

ieH  ie h 

X'  X' 

where  the  partial  thermal  conductivity  is  given  by  \'h  =  Mh{(f)hh ,  The  Erst 

term  of  the  heat  flux  represents  Fourrier’s  law,  the  second  term  corresponds  to  the  Dufour 
effect,  that  is,  heat  diffusion  due  to  the  partial  pressure  gradient,  electric  field,  and  average 
electron  force.  Finally,  the  third  term  the  transfer  of  heavy-particle  energy  due  to  diffu¬ 
sion.  In  Appendix  D,  we  show  that  the  heavy-particle  diffusion  velocities  and  heat  flux 
are  proportional  to  the  electron  driving  force  and  electron  temperature  gradient  through 
the  Fic  contribution  to  di,  i  G  H.  Kolesnikov  (1974)  has  already  introduced  electron 
heavy-particle  diffusion  coefficients  and  thermal  diffusion  coefficients  and  ratios  to  couple 
the  heavy-particle  diffusion  velocities  and  heat  flux  to  the  electron  forces.  Therefore,  we 
propose  to  refer  to  this  phenomenon  as  the  Kolesnikov  effect  for  the  heavy  particles.  Al¬ 
ternative  forms  of  the  heavy-particle  diffusion  velocities  and  heat  flux  are  also  introduced 

V,  =  -^^(dj  +  XjdJnT^,  ieH, 

ieH  .  ^  x  ^ 

qh  =  -A hdxTh  +  ph  xty  +  2^  PjhjVj , 

ieH  jeH 

where  the  thermal  diffusion  ratios  are  defined  from  the  relations 

E%y  =  y,  ien,  Eh  =  °. 

ieH  ieH 

and  the  thermal  conductivity  by  A h  =  A 'h  —  rth  Xj@j  ■ 

4.3  Second-order  electron  perturbation  function 

The  second-order  electron  perturbation  function  0e2  is  a  solution  to  the  electron  Boltzmann 
eq.  (22)  at  order  e°,  he., 

Te(0e2)  =  ^[-^(/e0,^)+3ee(/e0^,/e0^)+S' 

=  ~^:dxVh  -  SboPeV^-dl  - PeY l  ~ 

ieH 

where  4^,  /i  G  {rje,De,  H),  are  functions  of  Cc  (Graille  et  ah,  2009).  The  second- 
order  electron  diffusion  driving  force  is  defined  as  d2  =  —neqeM^VhAB/pe,  and  the  second- 
order  heavy-particle  driving  forces  d2  =  —  Vi:  i  G  II.  Uniqueness  of  the  solution  is  imposed 
by  means  of  the  constraints 

((/eVe2,^))e  =  0,  *  G  {1,2}, 


(59) 

(60) 
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i.e.,  the  perturbation  function  does  not  contribute  to  the  macroscopic  electron  mass  and 
energy.  By  linearity  and  isotropy  of  the  linearized  collision  operator  fFe,  the  second-order 
electron  perturbation  function  is  expressed  in  terms  of  the  heavy-particle  velocity  gradient, 
and  second-order  diffusion  driving  forces 

4>l  =  -4>le-dxvh  -  5mpc(j)^-dl  -  pe  (61) 

jeH 

where  the  vectorial  functions  0/e,  <p^e,  and  (f)^j  are  the  solutions  to  the  equations 
under  the  scalar  constraints 

((/e°^,^e))e  =  0,  J  G  {1,2}, 

with  n  G  {pe,  De,  (Dj)je H}- 


The  transport  fluxes  are  obtained  by  injecting  the  expression  for  02  given  in  eq.  (61) 
into  the  eqs.  (48),  (47),  and  (50).  We  obtain  the  second-order  electron  diffusion  velocity 
and  heat  flux 


Vc 2  =  -5b0Ded2e  -  J ^aejd 2,  q;  =  -5bopeded2e  -  pe  ^  Xjdj  +  PeKV2,  (62) 
je h  jeH 


and  the  average  electron  force  acting  on  i  heavy  particles 

Pe  ,  Pe 


F  = 


n,Mh 


<\,',dr 


riiMh 


XidxlnTe,  i  G  H. 


The  aei  coefficients  and  second-order  electron  thermal  diffusion  ratios  read 


a„  =  0»]„  xt  =  - ?,  •  6  H. 


(63) 


The  second-order  electron  diffusion  velocity  and  heat  flux  are  thus  proportional  to  the 
heavy-particle  diffusion  velocities,  that  is  the  Kolesnikov  effect  for  the  electrons.  To 
the  authors’s  knowledge,  it  is  the  first  time  that  such  second-order  transport  fluxes  are 
rigorously  derived  from  a  multiscale  analysis.  The  Kolesnikov  effect  terms  for  heavy 
particles  and  electrons  are  essential  to  derive  a  total  energy  equation  and  an  entropy 
equation  that  satisfy  respectively  the  first  and  second  laws  of  thermodynamics,  as  shown 
in  the  following  section.  Since  the  electron  collision  operator  is  of  the  order  of  1/e2  in 
the  electron  Boltzmann  eq.  (20),  it  is  important  to  mention  that  they  should  not  be 
confused  with  Burnett  transport  fluxes  (Ferziger  and  Kaper,  1972)  based  on  a  second- 
order  perturbation  function  and  a  collision  operator  of  the  order  of  1/e. 


5  Discussion 

5.1  Mass  conservation 

Summing  eq.  (41)  over  j  G  H  and  using  the  constraint  )Cjgh  Pj^j  =  6  given  in  eq.  (3.2), 
a  heavy-particle  mass  conservation  equation  is  obtained 

dtph  +  dx-(phvh)  =  0.  (64) 
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The  heavy-particle  mass  is  conserved  in  the  mean  heavy-particle  velocity  frame.  Then, 
adding  the  electron  drift  eq.  (39)  to  eq.  (64)  and  using  the  relation  pv  =  ptVh  +  £2pe[vh  + 
(V£  +  eVe2)/Mh],  a  conservation  equation  is  also  established 

dtp  +  dx-(pv)=  0,  (65) 

for  the  global  mass  p  =  ph  +  £2 pe-  We  notice  that  eq.  (65),  can  be  retrieved  by  summing 
the  Maxwell  transfer  eq.  (6)  over  the  species. 


5.2  First  law  of  thermodynamics 

A  flow  kinetic  energy  equation  is  obtained  by  projecting  eq.  (42)  onto  the  mean  heavy- 
particle  velocity 


dt(\ph\vh\2)  +  dx ■  vh{\ph\vh  |-  +  jj^p) 

L  1  h 

=  jtfiP  dx  vh  -  vh  dx  IJh  +  j^nqE-Vh  +  vh-Sb0I0AB. 


A  total  energy  equation  is  derived  by  adding  eqs.  (40)  and  (43)  to  the  flow  kinetic  energy 


dt(£)  +  dx •  (' Hvh )  =  -edx-{nh-vh)  -  A- dx  Q  +  I  E,  (66) 

where  quantity  £  =  pe  +  M%ph^\vh\2  stands  for  the  total  energy,  7i  =  £  +  p,  the  total 
enthalpy,  and  I  =  nqvb  +  [neqe(Ve  +  sV2)  +  £  JTeH  njQjVj\/Mh,  the  total  electric  current 
density.  The  term  IE  represents  the  power  developed  by  the  electromagnetic  field.  It 
has  the  form  prescribed  by  Poynting’s  theorem.  Hence,  the  first  law  of  thermodynamics 
is  satisfied.  We  notice  that  eq.  (66)  is  similar  to  the  Maxwell  transfer  eq.  (8)  for  the 
total  energy,  the  difference  is  attributed  to  the  heavy-particle  velocity  frame  chosen  for 
the  Chapman-Enskog  method. 


5.3  Second  law  of  thermodynamics 

In  addition  to  the  thermal  energy,  we  introduce  other  relevant  thermodynamic  functions. 
First,  the  species  Gibbs  free  energy  is  defined  by  the  relations 


pege  =  neTe  In 


nen 


o 


Te3/2Q° 


>  Pi9i  ^i-^h  111 


nm 


K:T,)3/2go 


where  the  translational  partition  functions  read 


/27rm°kBT°\3/2 

V  4  ) 


,  Q°h  = 


f  27rm?kBT0\  3/2 


Then,  the  species  enthalpy  is  given  by 


i  e  H, 


5  5 

Pe^e  -^neTe,  Pi^i  ^  ^  i  ^ 'h  ■  i  G  H. 
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Finally,  the  species  entropy  is  introduced  as 


s 


e 


i  G  H. 


Therefore,  the  mixture  entropy  reads  ps  =  YljeS  Pjsj-  The  thermodynamic  functions 
exhibit  a  wider  range  of  validity  than  in  classical  thermodynamics,  introduced  for  sta¬ 
tionary  homogeneous  equilibrium  states  (Giordano,  2008).  Indeed,  they  are  interpreted 
in  the  framework  of  kinetic  theory  by  establishing  a  relation  between  the  thermodynamic 
entropy  and  the  kinetic  entropy  (Giovangigli,  1999).  This  quantity  is  based  upon  the 
distribution  functions 


gkin 


/ 1  —  in 


'(27 r)3/2n°  ' 
rrijQl  3  _ 


dC,  + 


|l  —  In 


(27t)3'/277,° 

~qI  . 


d  cc. 


The  kinetic  entropy  and  the  thermodynamic  entropy  are  asymptotically  equal  at  order  e2 


Skin  =  ps  +  0{£2), 


provided  that  the  distribution  functions  follow  the  Enskog  expansion  given  in  eqs.  (22) 
and  (22).  Consequently,  a  first-order  conservation  equation  of  thermodynamic  entropy  can 
be  used  instead  of  a  conservation  equation  of  kinetic  entropy  to  ensure  that  the  second 
law  of  thermodynamics  is  satisfied.  A  global  entropy  equation  is  derived  in  Graille  et  al. 
(2009) 

dt(ps)  +  dx ■  (psvh)  +  dxJ  =  T,  (67) 

where  the  global  entropy  flux  is  given  by 


J  =  Jl  +  eJ\  +  eJl 


J<°  = 


MhTe 


Qe  Pe9eVe),  c  (<?e  Pe9e^)t  yJh  ~~ 


MhTe 


MhTh 


'<lh  -»^)- 

je  H 


For  weakly  magnetized  plasmas  (6  =  0)  and  unmagnetized  plasmas  ( b  <  0),  we  define 
Xh  =  ( dx\nTh,ph(di)ieu)T  and  xe  =  [dx\nTe,pe(de  +  £Sb0d2c)]T.  Hence,  the  global  entropy 
production  rate  reads 


( Te  —  T7)2  \  ^  Tij  1  .  .  1 

T  =  — — 2J  — +  £^S:S  +  £TT7r(AhXh'  ^  +  TTr 

TeTh  rrij  MhTh  MhTe 

jeH  J 


AeXei  Xg'j 


-£28m^Dedl-dl 


JbO 


Ml 


The  mass-energy  transport  matrices  Ah  and  Ae  dehned  in  Appendix  D  are  symmetric 
positive  semi-definite.  The  global  entropy  production  rate  is  thus  nonnegative  provided 
that  e  is  small  enough  in  the  6  =  0  case  and  that  the  collision  frequencies  uie,  i  G  H, 
are  nonnegative  in  the  two  6  =  0  and  6  <  0  cases;  the  second  law  of  thermodynamics 
is  satisfied.  The  Erst  term  of  the  global  entropy  production  rate  expresses  that  thermal 
nonequilibrium  between  the  electrons  and  heavy  particles  is  associated  with  a  production 
of  entropy.  When  Te  =  Th,  the  quasi-equilibrium  Maxwell-Boltzmann  distributions,  given 
in  eqs.  (25)-(27),  degenerate  to  the  equilibrium  Maxwell-Boltzmann  distribution,  given  in 
eq  (5),  prescribed  by  the  //-Theorem. 
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5.4  Onsager’s  reciprocal  relations 

In  this  section,  we  deduce  from  kinetic  theory  the  Onsager  reciprocal  relations.  The 
expressions  for  the  transport  fluxes,  denoted  by  the  vector  F,  are  proportional  to  the 
diffusion  forces,  denoted  by  the  vector  X,  i.e., 

Fa  ^  ^  L a/3  X^j. 

P 

Onsager’s  reciprocal  relations  are  symmetry  constraints  which  must  hold  between  the 
transport  coefficients  (see  de  Groot  and  Mazur,  1984;  Woods,  1986) 

La/3  =  [L^af-  (68) 

They  result  from  microscopic  reversibility.  We  identify  the  diffusion  forces  from  the 
quadratic  form  of  the  entropy  production  rate  given  in  eq.  (67)  and  use  the  transport 
coefficient  expressions  established  in  Section  4. 


At  order  £°,  the  first-order  electron  mass-energy  flux  vector 

F,  =  [9.  -  P.h.V„  VJT 

is  proportional  to  the  electron  diffusion  force  vector  Xc  =  (c^lnTe,  pede)T ,  as  expressed 
by  the  relation 

Fe  =  —  AeXc. 

The  matrix  Ae,  given  in  Appendix  D,  satisfies  Onsager’s  reciprocal  relations.  At  order 
el ,  the  momentum  flux  is  decoupled  from  the  mass  and  energy  fluxes  (de  Groot  and 
Mazur,  1984).  The  heavy-particle  viscous  tensor  flux  J7/(  is  proportional  to  the  force  S 
through  a  scalar  transport  coefficient,  the  shear  viscocity  Onsager’s  relation  is  trival. 
In  Appendix  D,  we  rewrite  the  mass  and  energy  transport  fluxes  in  terms  of  the  diffusion 
forces  by  replacing  expression  (63)  for  Fie  in  eq.  (54).  The  global  mass-energy  flux  vector 

F=  (ql-f>eKV?,  qk  —  PjhjVp  V2,  (VP)^ , 

ie  H 


is  proportional  to  the  global  diffusion  force  vector 

X  =  (dadnTe,  dx\nThl  phd'e ,  ph(d[)ie H)  , 
as  expressed  by  the  relation  F  =  —  A  X,  where  the  modified  driving  forces  read 


<  =  -d=,  < 

Ph 


—  dxPi  - 
Ph 


n,q,.E 

Ph 


i  G  H. 


Thus,  concerning  the  mass-energy  transport,  the  generalized  Onsager  reciprocal  relations 
for  the  second-order  electron  transport  coefficients  and  first-order  heavy-particle  transport 
coefficients  are  also  satisfied. 
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6  Conclusions 

In  the  present  study,  we  have  derived  from  kinetic  theory  a  unified  fluid  model  for  multi- 
component  plasmas  by  accounting  for  electromagnetic  field  influence,  neglecting  particle 
internal  energy  and  reactive  collisions.  This  model  is  valid  for  unmagnetized  and  weakly 
magnetized  plasmas  encountered  in  hypersonic  applications.  Given  the  strong  disparity 
of  mass  between  the  electrons  and  heavy  particles,  such  as  molecules,  atoms,  and  ions, 
we  have  conducted  a  dimensional  analysis  of  the  Boltzmann  equation,  following  Petit  and 
Darrozes  (1975),  and  introduced  a  scaling  based  on  the  £  parameter,  or  square  root  of  the 
ratio  of  the  electron  mass  to  a  characteristic  heavy-particle  mass.  The  multiscale  analysis 
occurs  at  three  levels:  in  the  kinetic  equations,  the  collisional  invariants,  and  the  collision 
operators.  The  Boltzmann  equation  has  been  expressed  in  the  heavy-particle  reference 
frame  allowing  for  the  first-  and  second-order  electron  perturbation  function  equations 
to  be  solved,  as  opposed  the  inertial  reference  frame  chosen  by  Degond  and  Lucquin- 
Desreux  (1996a, b).  The  system  has  been  examined  at  successive  orders  of  approximation 
by  means  of  a  generalized  Chapman-Enskog  method.  The  micro-  and  macroscopic  equa¬ 
tions  derived  at  each  order  are  reviewed  in  Table  2.  Depending  on  the  type  of  species, 
the  quasi-equilibrium  solutions  are  Maxwell-Boltzmann  velocity  distribution  functions  at 
either  the  electron  temperature  or  the  heavy-particle  temperature,  thereby,  allowing  for 
thermal  nonequilibrium  to  occur.  At  order  e1,  the  set  of  macroscopic  conservation  equa¬ 
tions  of  mass,  momentum,  and  energy  comprises  multicomponent  Navier-Stokes  equations 
for  the  heavy  particles  and  first-order  drift-diffusion  equations  for  the  electrons.  The  ex¬ 
pressions  for  the  transport  fluxes  have  also  been  derived:  first-  and  second-order  diffusion 
velocity  and  heat  flux  for  the  electrons,  and  first-order  diffusion  velocities,  heat  flux, 
and  viscous  tensor  for  the  heavy  particles.  The  transport  coefficients  have  been  writ¬ 
ten  in  terms  of  bracket  operators.  We  have  also  proposed  a  complete  description  of  the 
Kolesnikov  effect,  i.e.,  the  crossed  contributions  to  the  mass  and  energy  transport  fluxes 
coupling  the  electrons  and  heavy  particles.  This  effect,  appearing  in  multicomponent 
plasmas,  is  essential  to  obtain  a  positive  entropy  production.  It  also  contains,  as  a  de¬ 
generate  case,  the  single  heavy-species  plasmas  considered  by  Degond  and  Lucquin  for 
which  the  Kolesnikov  effect  is  not  present.  The  properties  of  symmetry  and  positivity 
of  the  electron  and  heavy-particle  mass-energy  transport  matrices  imply  that  the  second 
law  of  thermodynamics  is  satisfied,  as  shown  by  deriving  an  entropy  equation.  Moreover, 
Onsager’s  reciprocal  relations  hold  between  the  transport  coefficients.  The  first  law  of 
thermodynamics  was  also  verified  by  deriving  a  total  energy  equation.  Finally,  the  purely 
convective  system  of  equations  is  found  to  be  hyperbolic  in  Graillc  et  al.  (2009),  thus 
leading  to  a  well  defined  structure.  The  explicit  expression  for  the  diffusion  coefficients, 
thermal  diffusion  coefficients,  viscosity,  and  partial  thermal  conductivities  can  be  obtained 
by  means  of  a  variational  procedure  to  solve  the  integral  equations  (Galerkin  spectral 
method,  see  Chapman  and  Cowling,  1939).  The  expressions  for  the  thermal  conductivity, 
thermal  diffusion  ratios,  and  Stefan-Maxwell  equations  for  the  diffusion  velocities  can  be 
derived  by  means  of  a  Goldstein  expansion  of  the  perturbation  function,  as  proposed  by 
Kolesnikov  and  Tirskiy  (1984).  Finally,  the  mathematical  structure  of  the  transport  ma¬ 
trices  obtained  by  the  variational  procedure  can  readily  be  used  to  build  efficient  transport 
algorithms,  as  already  shown  by  Ern  and  Giovangigli  (1994)  for  neutral  gases,  and  Magin 
and  Degrez  (2004a)  for  unmagnetized  plasmas. 
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A  H-Theorem 

We  consider  here  the  simplified  case  of  a  spatially  uniform  gas  in  the  absence  of  external 
forces.  In  this  case,  the  Boltzmann  eq.  (2)  reads 

a.*/*  =  E3«h*-h)-  ieS- 

jeS 

By  taking  the  time  derivative  of  this  quantity  and  using  the  simplified  Boltzmann  equation 

^  =  £  /(In/Z  +  DA-Vdc* 

ies  J 

=  £  /  (/.*T -/.*//)  (In /.*  +  l)K-yiydu,dc*dc* 

ties  J 

=  £  [  (/.*£' -/.*//)  (Ill //  +  l)|c*-c*|<r*  du,dc*dc* 

ties  J 

=  £  [  (/.*//- //*//*)  (In//*  +  l)|c'*-c“|ff*jdWdc'*dc;* 
ties  J 

=  £  [ (ffff  -//*/;*)  (iny  +  i)|c'*-c;*i^da,dc'*dc;* 

ties  J 

=  -j£  /  U*7/'-ii*/;)[ln(//*/;*)-ln(/.*//)]|c*-c*|ff*,dWdc*dc* 
ties  J 

where  we  have  used  the  equalities  dc'*dc'*  =  dc*dc*  and  c'*  —  c'*|  =  |c*  —  c*|. 


B  Electron  heavy- particle  interactions 


The  study  of  the  electron  heavy-particle  collision  dynamics  yields  the  dependence  of  the 
peculiar  velocities  on  the  e  parameter.  First,  we  express  momentum  conservation  in  terms 
of  the  peculiar  velocities  in  the  heavy-particle  velocity  frame.  Considering  a  collision  of  a 
heavy  species,  «  6  H,  against  an  electron,  the  peculiar  velocities  after  collision  C[  and  C'e 
are  related  to  their  counterpart  before  collision  C,;  and  Cc 


C[  =  - +  mi  2  C‘  +S  - 

rrii  +  el  rrii  +  e2  rrii  +  £2 


£ 


fi'  - 

v-ye  i  o  i  .  o 

m*  +  77li  +  £2 


cc  +  ^^a-s  m 


rrii  +  £2 


|e  Q  —  Ce  u>,  i  6  H, 

|  £  —  C*e  |  CJ , 


provided  that  the  mean  heavy-particle  velocity  is  not  modified  by  this  single  collision 
event.  The  direction  of  the  relative  velocities  after  collision  is  defined  in  their  center  of 
mass  by 


uj  = 


sC'-C' 
eC[  -  C' 
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Symbol  s  stands  for  an  integer  either  equal  to  +1  for  the  collision  operator  die,  i  E  H,  or 
—  1  for  Sei,  i  E  H.  The  crossed-collision  operators  are  expanded  in  Graille  et  al.  (2009). 
The  collision  operator  die,  i  £  H,  can  be  expanded  in  the  form 


&U./e)(a)  =  +e23Ufi,fe)(C,)  +  e33Ufi,f,)(Ci)  +  0(e*), 


where  the  zero-order  collision  operator  d%(fi ,  fe)(Ci),  «  6  H,  vanishes.  The  collision 
operator  3ei,  i  G  H,  can  be  expanded  in  the  form 


a„(/,./i)(a)  =  3i(/„/i)(c'.)+£ai(/,./,)(a)+e23i(/../i)(c,)+e3ai(/.,/i)(c'«)+c’(e4). 


We  also  define  a  collision  frequency  as  an  Maxwell- Boltzmann  averaged  momentum 
cross-section 

i '<•  =  ^  / Qg’tic.rjicyftc.)  dc„,  i  €  h, 

and  introduce  the  generalized  momentum  cross-section  of  Chapman  and  Cowling  (1939) 
in  a  thermal  nonequilibrium  context  (Mitchner  and  Kruger,  1973) 

/»7r 

Q[l\\Ce\2)  —  2n  /  o-,:e(|Ce|2,  cos6*)(l  —  cos#)  sin#  d#,  i  6  H, 

Jo 

where  symbol  9  stands  for  the  angle  between  the  vectors  uj  and  e. 

C  Whale  equations 

Based  upon  the  dimensional  analysis  of  Section  3.2,  the  electron  Boltzmann  eq.  (20) 
becomes 

£“1^e“1(/e°>  <k)  +  @e(fe,  4,  <Pe)+£^e(fe,  4, 4>l,  <t>e)  =  ^d^+^d^+de+Edl+O^2), 

where  the  electron  streaming  operators  read  at  successive  orders 


«,-*(. /«",*)  =  ifca-o*/,"  +  qc(EtE + sKc^B)dcj,°, 

»“(/<“, 4, 42)  =  a,f°  +  ^c.-s.u'Vb)  +  vh-dxf®  -  (ocJ,°®c,):dxvh 

+  9e{Pb0^hvh^B  +  5b(^i)CeABydcBfc  +  +  ^boCe^B)  ■  dce  ) ) 

9j(/,°,4,tf,ri?)  =  dt(fe%)  +  jECe.dx(fX)  +  vh-dx(fe\)  -  Mb^-dcJ 

-  ( dcB{f,J(jk)®Cc)-.dxvh  +  qe(Sb(-i)MhvhAB  +  5b{_2)CeAB) -dcJe0 

+  qe(Sb0MhvhAB  +  5b(-i)CeAB)-dce(f%)  +  ^(AE  +  Sb0CeAB)-dcM 

The  electron  collision  operators  are  given  by 

3;2 

je  H 

s;'  =  3«(/,°4  ,f°)+  3«(/=°,  f°4 k )  +  E  4”)  +  + iiyfrfl, 

ieH 


eVe) 
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t  =  3==(/eVe  ,  /<")  +  3Uf,% ,  f°A )  +  3 «(/«",  /e°42)  +  E  3°,  (/„V„2,  //)  +  +  3! 


ieH 


3°  =  E  .  li*i  )  +  +  dljife,  //*  )  +  3i(/«°,  //), 

ieH 

^  =  aee(/eVe3,  fe)  +  3ee(fUe,  A  )  +  3ee(fe%  ,  /eVe)  +  See(fe°,  fUe) 

+  E  ^(/eVe,  ff)  +  +  £, 


ie  H 


3;  =  E  { H )  +  3°i(A°4,//#)  +  +  3p/e°4,//4) 

+  3p/„°,  f°$)  +  3  i(/„"4 ,  //)  +  3i(/e°,  )  +  3i(/„°,  4°)} . 


ieH 


For  ease  of  readability,  we  strike  through  the  collision  operators  that  vanish  when  /e°  and 
4  !  6  H,  are  isotropic  functions.  Likewise,  the  heavy-particle  Boltzmann  eq.  (21)  is 
found  to  be 


4(f)  +  e^lUiA)  =  s-%1  +  a°  +  4  +  0(e2),  i  E  H, 
where  the  heavy-particle  streaming  operators  read  at  successive  orders 

% Hf°A)  =  W%)  +  {-irhC,+vh)-dM%)  +  Sfr>[(Ci+Mhvh)fJi].9ctj? 

+  ^kE-9c,(m-  M^-ea(j%)~  (9c i(m^Ci)-9.vh. 

The  heavy-particle  collision  operators  are  given  by 

3r'  = 


jeH 


3“  =  E  3«(/i°4 .  //)  +  3«(/i°,  f°4>j )  +  +  3”, 

ieH 

a°  =  aff /eVe)  +  al(f  f), 

a!  =  E  Mf  f ,  f )  +  (f  * » f:% )  +  Mf  f  f )  +  alif^ffT + a), 


J'GH 


a)  =  aff  ,  /eVe )  +  aL(f  ff)  +  sUf% ,  f)  +  af  f ,  fe% )  +ii44T- 


D  Transport  matrices 

The  first-order  heavy-particle  mass-energy  flux  vector  Fh  =  (qh  —  X)jeH  PjhjVji  CK)» sh)T, 
is  proportional  to  the  heavy-particle  diffusion  force  vector  Xh  =  ( dxlnTh ,  Ph(di)ieE)T ,  as 
expressed  by  the  relation  Fh  =  —AhXh,  where  the  heavy-particle  mass-energy  transport 
matrix  is  given  by 

,  _  (  ThXh  [(9f)teHF 
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The  first-order  electron  mass-energy  flux  vector  Fe  =  [qe  —  peheVe,  V^T  is  proportional 
to  the  electron  diffusion  force  vector  Xc  =  (c^lnTe,  pede)T,  as  expressed  by  the  relation 
Fe  =  —  AeXc,  where  the  electron  mass-energy  transport  matrix  is  given  by 


Ac  = 


(TeK 

V 


The  heavy-particle  diffusion  velocities  given  in  eq.  (56)  read 


Vt  =  -0?djn  Te 


e?dJnTh 


Diede 


E  Wi • 

j'e  H 


GH, 


with  the  modified  driving  forces 


7/  Pe  i  i/  1  rj  niQi  j-,  t  t 

dc  =  —de,  d-  =  — Oxpi - E,  i  G  H. 


Ph 


Ph 


Ph 


The  matrices  of  heavy-particle  electron  diffusion  coefficients  and  electron  thermal  diffusion 
coefficients  are  defined  as 

a.  =  XX«*.  =  ieH- 

iSH  Ph  jgH 

The  heavy-particle  heat  flux  given  in  eq.  (58)  reads 

qh  =  -A  hc'dxTe  -  \’hdxTh  -  ph6hed'c  -phJ2  +  Y  PihiVv 

jSH  jSH 

with  the  matrices  of  partial  thermal  conductivity  and  thermal  diffusion  coefficient 


JjOley 


ieH 


jGH 


Then,  substituting  the  previous  expression  for  \P,  i  G  H,  into  the  eq.  (62)  for  the  second- 
order  electron  diffusion  velocity  Vc2,  one  has 

=  -6Pedx\nTe  -  9hcdxlnTh  -  Deed'e  -  Y  De]d'p 

tGH 

with  the  following  second-order  matrices  of  electron  diffusion  coefficients,  electron  heavy- 
particle  diffusion  coefficients,  and  electron  thermal  diffusion  coefficient 

Dee  =  Y  ®ejDje,  Dei  =  Die,  i  G  II,  0ee  =  Y  aejOej. 

je H  ieH 

Regarding  the  electron  heat  flux  given  in  eq.  (62),  one  obtains 


ql  =  -A  l'dxTe  -  pA  hJdxTh  -  ph9eed'e  ~PhY  °eA  +  PeheVe2’ 

h  jeH 
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with  the  second-order  matrix  of  electron  partial  thermal  conductivity  Xf  =  ne  ^  -gH  XjQj- 
The  global  mass-energy  flux  vector  F  =  (q2e- peheVe2 ,  qh~Y,j&n  PjhjVji  Ve 2,  (K)*eH)T,  is 
proportional  to  the  global  diffusion  force  vector  X  =  (dx\nTei  dxlnTh,  Ph,d'e,  h)T, 

as  expressed  by  the  relation  F  =  —A  X,  where  the  mass-energy  transport  matrix  has  the 
following  block  structure 


/(A?re)  (a yzy 
(A  JT,)  {\'„Th) 


H  H 
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